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INTRODUCTION 

A substantial  amount  of  effort  has  gone  into  predicting  the  added  mass  and 
damping  coefficients  of  an  oscillating  ship.  Historically,  a large  part  of 
this  effort  has  been  concentrated  on  the  vertical  plane  motions  of  heave  and 
pitch  with  some  interest  shown  in  the  horizontal  motion  of  roll.  Korvin- 
Kroukovsksy  and  Jacobs  (1957)  emphasized  the  importance  of  coupling  coeffi- 
cients in  the  vertical  plane  motions,  i.e.  the  pitch  force  due  a heave  motion 
and  vice  versa.  Their  method  of  calculating  these  coefficients  was  based  on 
strip  theory  with  some  adjustments  made  to  include  forward  speed  and  three- 
dimensional  effects.  However,  their  method  for  finding  these  "dynamic  coupling" 
terms  appears  to  be  incomplete  in  that  these  terms  do  not  satisfy  the  symmetry 
relations  established  by  Timman  and  Newman  (1962) . 

Salvesen,  Tuck  and  Faltinsen  (1970)  applied  a more  consistant  approach  to 
the  problem  of  ship  motions  and  developed  a theory  for  both  the  vertical  and 
horizontal  modes  of  motion.  Their  cross-coupling  coefficients  did  satisfy  the 
Timman-Newman  (1962)  relations  and  gave  improved  results  when  theory  was  com- 
pared with  experiment. 

Ogilvie  and  Tuck  (1969)  found  the  added  mass  and  damping  coefficients  for 
heave  and  pitch  by  using  a systematic  applic  ation  of  matched  asymptotic  expan- 
sions. Their  results  satisfied  the  Timman-Newman  (1962)  symmetry  relations 
but  differed  from  the  coefficients  given  by  Salvesen,  Tuck  and  Faltinsen  (1970) 
in  a number  of  ways.  The  Ogilvie-Tuck  coefficients  included  a term  which  rep- 
resented the  integral  of  the  square  of  the  velocity  potential  evaluated  on  the 
free  surface.  They  did  not  include  a velocity-squared  term  which  they  consid- 
ered to  be  of  higher  order.  In  an  effort  to  determine  the  relative  importance 
of  the  free  surface  integral  terms,  Faltinsen  (1974)  evaluated  them  and  compared 
them  with  both  experiments  and  the  terms  given  by  Salvesen,  Tuck  and  Faltinsen. 
The  results  indicate  that  the  Ogilvie-Tuck  heave-pitch  coupling  coefficients 
are  important  and  compare  better  with  experiment  than  previous  theories. 

Timman  and  Newman  (1962)  included  the  horizontal  motions  in  their  sym- 
metry relations.  Specfically  they  reported  that  the  cross-coupling  terms 
between  roll  and  yaw  and  yaw  and  sway  were  antisymmetric  with  respect  to  for- 
ward velocity.  Inspection  of  the  Salvesen,  Tuck  and  Faltinsen  coefficients 
show  that  they  satisfy  these  conditions. 
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One  motion  coefficient  in  the  horizontal  plane  that  has  received  some 
attention  is  the  roll  damping  coefficient.  Typically  the  damping  coeffi- 
cient is  composed  of  a velocity  independent  part,  calculated  from  potential 
theory,  and  a correction  factor  used  to  account  for  viscous  effects.  In  an 
effort  to  more  clearly  understand  these  various  effects  on  roll  motion,  Sugai 
and  Yamanouchi  (1963)  conducted  a series  of  experiments  using  a self-propelled 
model  with  opposing  gyroscopes  to  provide  a rolling  moment.  One  of  the  re- 
sults of  these  experiments  was  the  indication  that  rolling  motion  becomes 
more  linear  as  forward  speed  is  increased.  In  other  words,  the  non-dimen- 
sional quantity  of  roll  motion  divided  by  roll  exciting  moment  (multiplied 
by  the  appropriate  constant  to  correct  for  the  units)  becomes  less  a function 
of  the  actual  magnitude  of  the  rolling  moment.  This  implies  that  roll  damping 
is  speed  dependent  and  that  this  dependency  may  be  at  least  as  important  as 
the  viscous  damping.  Watanabe  (1977)  applied  the  principle  of  thin  ship 
theory  to  this  problem  and  produced  results  that  show  a speed  dependency  on 
the  roll  damping  coefficient.  However,  he  did  not  conclusively  state  what 
effects  a non-thin  ship  would  have  on  his  theory. 

Using  the  same  assumptions  made  by  Ogilvie  and  Tuck  (1969)  for  the  ver- 
tical plane  motions,  we  will  consider  in  this  paper  the  following  two  ques- 
tions: First,  will  a consistant  slender  body  theory  produce  roll  damping 

coefficients  that  reduce  roll  motions  as  the  speed  of  the  ship  increases? 

And  second,  will  the  theory  produce  numerically  significant  terms  to  the  ones 
derived  by  Ogilvie  and  Tuck. 


■ 


STATEMENT  OF  THE  PROBLEM  WITH  RESULTS 


equations  OF  MOTION 

Similar  to  that  shown  by  Salvesen,  Tuck  and  Faltinsen  (1970) , the  coupled 
equations  of  motions  for  sway,  roll,  and  yaw  may  be  given  as  follows: 

(A22+  M)H2  + *22*2  + (A24-  MZC)  n4  + B24% 

+ (A26+  XcM)H6  + B26\  = U) 

(A42-  Mzc>  V B42V  (A44+  VH4  + B44^4  + C44n4 

+ ( A - I ) n,  + B n = f ela)t  (2) 

46  46  6 46  6 4 

ia62+  xc„)i;2  ♦ b62j2  ♦ <a64-  i46)  n4  ♦ b64;4 

+ (A66+  I6)^6  + B66\  = F6e  (3> 

where 

M is  the  mass  of  the  ship 

Ajk'Bjk  are  mass  and  damping  coefficients  respectively 

I . is  the  moment  of  inertia  about  the  j-th  axis 
3 J 

Ijk  is  the  product  of  inertia 

is  the  vertical  location  of  the  center  of  gravity 
iwt 

F_.e  is  the  force  or  moment  in  the  j-th  mode  due  to  waves 
C44  is  the  hydrostatic  roll  restoring  moment 

*c  is  the  longitudinal  location  of  the  center  of  gravity 

oj  is  the  frequency  of  encounter 

for  j=2 ,4,6  is  the  sway,  roll,  and  yaw  displacement,  respectively. 

(The  dots  denote  time  derivatives,  i.e.  n.  is  the  roll  accel- 
eration.) 

The  coordinate  system  used  is  a right  hand  one  with  the  origin  located  in 
the  plane  of  the  undisturbed  free  surface  and  the  z axis  passing  through  mid- 
ship. See  Figure  1 for  definitions  of  positive  sway,  roll,  and  yaw. 
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Figure  1 Coordinate  System 


Consider  the  coordinate  system  as  fixed  in  an  incident  stream  with  veloc- 
ity U flowing  in  the  positive  x direction  and  the  ship  swaying,  rolling, 
and  yawing  about  that  system.  We  now  desire  to  find  the  added  mass  and  damping 


coefficients,  A., 
Ik 


and  B 


jk 


for  sway,  roll,  and  yaw  using  assumptions  sim- 


ilar to  the  ones  used  by  Ogilvie  and  Tuck  (1969) . 


If  we  define  the  motion  of  the  ship  in  a two  parameter  expansion,  e , 
a slenderness  parameter  related  to  the  beam  to  length  ratio,  and  6 , a 

motion-amplitude  parameter  related  to  the  smallness  of  motion,  we  can  require 
that  the  motion  be  smaller  than  the  beam  of  the  ship  even  as  e->0  . Specifi- 

cally, assume  that 

displacements  resulting  from  n.  = 0(6e) 


where  is  the  motion  in  the  j-th  mode.  Also  assume  that  the  frequency  of 

encounter  is  of  the  following  order: 


0(C-1/2) 


Then  the  velocity  potential  for  the  complete  (linearized)  solution  can  be 
represented  as  follows: 

<)>(x,y,z,t)  = Ux  + UX(x,y,z)  + iMx,y,z,t)  (4) 


where  the  first  two  terms  give  the  solution  of  the  steady-motion  problem  as 
shown  by  Tuck  (1965)  and  the  last  term  represent  everything  that  must  be  ad- 
ded in  order  to  satisfy  the  boundary  conditions.  We  assume  that  iHx,y,z,t) 
has  a time  dependence  of  elwt  and  the  velocity  U is  of  order  one. 

As  shown  by  Ogilvie  and  Tuck  (1969) , we  may  put  equation  (4)  into  the 
governing  equations  describing  the  boundary  value  problem  of  the  oscillating 
ship.  We  next  linearize  the  problem  with  respect  to  the  amplitude  motion,  but 
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keep  higher  order  terms  of  the  slenderness  parameter  expansion.  The  time-depen- 
dent part  of  the  potential  then  is  found  to  satisfy  the  following  equations: 

ip  + ip  = 0 in  the  fluid  domain,  (5) 

yy  zz 

—(i)2  ip  + gip  = -ia)U(2iJ)  + 2y  ip  + X <l>)  on  z=0,  (6) 

z x y y yy 

and  on  the  hull 


3^  = iwn2(n2+  xn&)  + i<i)n4(yn3-  zn2)  + un2n6 

-U(n2+  ^n6)(n2xyy+  n3xyz)  + up4[(n3xy-  n2xz) 

+ z(n2Xyy+  Vyz*  " y(n2Xyz+  n3Xzz^  (7) 

where  n is  the  unit  normal  directed  out  of  the  fluid,  n2  and  n3  are  the 
components  of  the  unit  normal  in  the  y and  z direction  respectively  and 
varibles  subscripted  with  y and/or  z denote  partial  differentiation  with 
respect  to  that  coordinate. 

The  form  of  the  above  equations  (5)  - (7)  can  only  be  correct  if  the  for- 
mal rules  of  matched  asymptotic  expansions  are  followed.  In  other  words,  the 

far  field  expansion  of  the  <|<(x,y,z,t)  potential  must  match,  to  an  appropriate 
order,  an  inner  expansion  of  a potential  representing  a line  of  pulsating 
sources  and  dipoles.  To  show  that  this  is  indeed  the  case,  we  could  use  a 
method  similar  to  that  used  by  Troesch  (1975)  , where  he  applied  the  theory 
of  Fourier  Transforms  or  use  the  method  of  complex  variables  as  shown  by 
Ogilvie  (1974).  In  both  cases,  the  authors  were  solving  anti  symmetric  prob- 
lems, which  are  applicable  to  the  sway,  roll,  and  yaw  potential  discussed  in 
this  paper. 

A solution  for  the  complete  i|)(x,y,z,t)  problem  can  now  be  given  in  much 
the  same  manner  as  shown  by  Troesch  (1975)  for  the  case  of  an  anti  symmetric 

pressure  distribution  on  the  free  surface.  The  actual  solution  is  fairly  com- 

plex and  it  is  not  necessary  to  repeat  it  here.  Rather,  we  note  that  since  a 
solution  exists,  we  may  find  the  pressure  from  Bernoulli's  equation  and  sub- 
sequently the  hydrodynamic  force  acting  on  an  oscillating  ship.  The  details 
are  given  in  the  appendix. 

ADDED  MASS  AND  DAMPING  COEFFICIENTS 

Following  Ogilvie  and  Tuck  (1969) , the  generalized  hydrodynamic  force, 


F j , may  be  related  to  the  added  mass  and  damping  coefficients  in  the  fol- 
lowing manner 

Fj(t>  - £&2v  “‘jkHV'15  * Tji’  * Tjk’  >k(t' 


where  the  superscript  (0)  denotes  the  usual  zero  speed  strip  theory  terms. 

The  velocity  dependent  terms,  which  are  of  higher  order  are  found  in  T^1)  and 
. These  coefficients  can  then  be  given  in  terms  of  the  steady  motion 
potential,  x » and  the  usual  zero  speed  forced  oscillation  potential,  ^ , 

defined  in  the  following  manner: 

<J> . + <f> . =0  in  the  fluid  domain 

j 1 

yy  zz 

9<f> . 

~ N.  > j = 2,4,6,  on  the  hull 
oN  ] 

-ai2ij) .+  g<p.  = 0 on  z = 0 
3z 


and 


!>  .-*■  + A .e 
1-3 


riv \- 


as  y •+  + oo 


where  N is  now  the  two  dimensional  normal  in  the  y-z  plane  and  N_.  is 
given  as  follows: 

N2  " n2 


N4  = n3y 


xn. 


V 


The  complete  expressions  for  the  complex  factors  T_.^  are  given  in  the 
Appendix.  If  we  make  two  assumptions,  we  can  simplify  the  expressions  for 
T24^  ' T24^  ' T46^  ' and  considerably.  First  assume  that  the  ship  is 


46 


symmetrical.  (For  a symmetrical  ship  the  half-beam  at  the  waterline,  yQ(x)  • 
is  an  even  function  with  respect  to  x and  the  y and  z derivatives  of  the 
steady  motion  potential,  x(*/Y»z)  , are  odd  functions  with  respect  to  x .) 

Next  assume  that  the  ship  in  question  has  long  sections  of  parallel  mid-body, 
that  is,  long  relative  to  its  end  section.  (This  has  the  effect  of  making  the 
end  contributions  to  T and  T j2)  higher  order  than  the  mid-body  contribu- 
tions.) The  complete  factors  T may  then  be  written  in  the  following  manner: 

jK 


] 


A 


r 

1 


and 


T22^  = “P(iu)2 


dx 
L C (x) 


d*n2*2  = 


dx  t 


(0) 

22 


= -p  ( ito)  2 | dx  | d«.(n3y  - n2z)  ^ 

T<»l  - { dx  X* 

L 

T<J}  = T<°>  = -p(ico)2  j dx  |dln2^ 


L C (x) 


,2  t(0) 

22 


T26}  = T62  - { dx  x 
L 


T46>  = T64>  = -p(iw): 


= 0 , j = 2,4,6 
DD 

TO)  = -t(')  = — T(o) 
26  62  iw  22 

T(l)  = TU)  = o 

T(i)  _ T(i)=  JL  T(o) 

64  ~ 46  ico  42 

T<2)  = 0 , j = 2,4,6 

n 

t(2)  = T (2 ) = 0 
24  42 


L C(x) 


t(°) 

Z22 


dx  x 


d*V4 


C (x) 


r26J  = “T62)  = p{iu,)3  ^ j d*[  j dy(^2  " A2e_2iVy) 

L yo 

i A2  e-2ivy0] 


T(2) 
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= = p(iw)3  T’  j <**[  j dy(^2V  A2A4e' 

L yo 

- -h  V4*'21vy°] 


■2ivy 


r — 


8 


where 


p is  the  water  density 

C(x)  is  the  hull  contour  at  station  x 

L is  the  ship  length 

v is  the  wave  number  found  from  ui2  = vg 
g is  the  gravitational  constant 
y0 (x)  is  the  half-beam  at  station  x 


—ivy 

A.e  is  the  behavior  of  d> . as  y-x». 

1 1 

The  added  mass  and  damping  coefficients  are  then  given  by  the  relations 

A.,  = \ Re  {t(°>  + t!.1)  + T<2)  } (8) 

jk  a)'1  3k  3K  3k 


B.,  = — Jm  + t(3)  + T<2>)  (9) 

3k  u>  3k  3k  3k 

where  He  and  Im  denote  the  real  and  imaginary  parts  of  the  complex  expres- 
sion respectively. 

Before  we  compare  these  coefficients  with  those  derived  by  Salvesen,  Tuck 

and  Faltinsen  (1970),  write  Tj(2)  and  t{2^  following  a notation  established 

ib  46 

by  Ogilvie  and  Tuck  (1969)  and  used  by  Faltinsen  (1974) . Let 
P ( iw)  3 ^ j ds  S’2  = T<2) 


p(ia))3  j ds  *2*4  E T46> 


where  a bar  has  been  drawn  through  the  integral  sign  to  indicate  that  the  in- 
tegral does  not  really  exist  as  written,  and  F denotes  that  the  limits  of 
integration  extend  on  the  free  surface  from  the  body  to  infinity. 


The  added  mass  and  damping  coefficients  derived  by  the  two  different  methods 
are  given  in  Table  1.  The  ship  is  assumed  to  be  pointed  at  both  ends  so  that 
the  end  terms  given  in  the  Salvesen*  Tuck,  and  Faltinsen  (1970)  coefficients 
do  not  appear.  Also  a^  and  b^  represent  the  usual  two  dimensional  sec- 
tional added  mass  and  damping  coefficients. 
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NUMERICAL  RESULTS 


Recall  that  we  are  addressing  ourselves  to  two  questions.  First,  will 
this  theory  show  a n|*>o<i  dr>|*>ii(l<*nc<'  for  the  roll  damping  coefficient  and,  sec- 
ond, are  the  extra  terms,  Tjk^  * numerically  significant?  In  order  to  answer 
the  first  question,  set  F_  in  equation  (1)  and  F,  in  equation  (3)  equal 

2 b 

to  zero  and  let  F4  in  equation  (2)  be  some  constant.  This  corresponds  to 

the  experiments  described  by  Sugai  and  Yamanouchi  (1963)  where  they  applied 

a pure  roll  moment  to  a moving  model.  A computer  program  was  written  that 

solved  equation  (1)  - (3)  using  the  added  mass  and  damping  coefficients  given 

in  equations  (8)  and  (9) . The  hull  offsets  used  were  from  a Series  60, 

C = .70  hull  and  the  model  was  assumed  to  have  the  following  additional  char- 
B 

acteristics : 

z = -0. 2d 
c 

I . = Ml~(0.397B)2  + (z  ) 2 ~1 
4 “ C 

C44  = MgGM 
GM  = -z 

C 

I,  = M(0.25L)2 
6 

I46  = °‘° 

where  d is  the  draft  of  the  ship,  B is  the  beam,  and  L is  the  length. 

The  computer  program  that  calculated  the  horizontal  plane  coefficients 
also  returned  the  ones  for  the  vertical  plane,  i.e.,  the  Ogilvie-Tuck  added 
mass  and  damping  coefficients  for  heave  and  pitch.  The  Series  60  hull  form 
was  selected  in  order  to  compare  the  results  returned  for  heave  and  pitch 
with  those  reported  by  Faltinsen  (1974).  (The  comparison  was  satisfactory.) 
While  the  hull  form  selected  is  not  truly  symmetrical  and  does  not  have  large 
sections  of  parallel  mid-body,  it  should  satisfy  the  conditions  of  symmetry 
and  relatively  constant  cross  sections  to  a sufficient  degree  to  answer  the 
two  questions  posed  in  the  Introduction. 

The  method  used  for  evaluating  the  free  surface  integrals  is  described 

in  the  Appendix.  The  results  of  the  forced  rolling  tests  for  three  different 

Froude  numbers  (F  = .15,. 20,. 30)  are  shown  in  Figures  2,3,  and  4.  Inspec- 
n 

tion  of  Table  1 shows  that  b44  is  independent  of  speed.  Therefore,  the  only 
way  that  the  roll  motion  at  resonance  will  be  influenced  by  speed  is  through 
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the  coupling  between  roll  and  sway  and  roll  and  yaw.  Figure  3 indicates  that 
this  effect  is  negligible.  Only  the  motion  of  yaw,  in  Figure  4,  shows  any 
speed  effects.  We  can  conclude  that  slender  body  theory,  as  derived  in  this 
paper,  will  not  show  an  increase  in  the  roll  damping  coefficient  as  forward 
speed  increases.  This  result  can  be  contrasted  with  the  results  given  by 
Watanabe  (1977)  where  an  application  of  thin-ship  theory  did  produce  larger 
roll  damping  coefficients  as  the  forward  speed  increased.  It  should  be  noted, 
however,  that  the  theory  of  that  paper  included  a trailing  vortex  sheet  and 
consequently  differed  fundamentally  from  the  theory  presented  here.  In  a 
recent  paper  by  Ikeda,  Himeno  and  Tanaka  (1978) , it  is  suggested  that  the 
"wave  damping  component"  of  a moving,  rolling  ship  is,  to  a large  extent. 


independent  of  speed  for  values  of  t greater  than  0.5.  Here  t equals 
Ucj/g  . The  "wave  damping  component"  corresponds  to  the  damping  coefficient 
B44  presented  in  this  paper.  Our  assumptions  make  t=0(e-1/2)  , which 

represents  high  speed  and/or  high  frequency  of  rolling.  Thus  it  appears  that 
is  consistant  with  the  results  presented  by  Ikeda,  et.al.  (1978). 

To  demonstrate  the  difference  between  the  added  mass  and  damping  coef- 
ficients as  derived  by  the  theory  presented  in  this  paper  and  those  coeffi- 
cients derived  by  Salvesen,  Tuck  and  Faltinsen  (1970) , Figures  5 through  8 


are  presented.  The  coupled  sway-yaw  added  mass  and  damping  coefficients  are 

non-dimensionalized  by  pVL  and  pVL/g/L  respectively.  Here  V is  the 

volume  displacement  of  the  ship.  The  coupled  roll-yaw  added  mass  and  damping 

coefficients  are  non-dimensionalized  by  pVL2  and  pVL2 /g/L  respectively. 

They  are  all  plotted  as  functions  of  w/L/g  . The  results  are  shown  for 

a F =0.2  only.  Froude  numbers  of  0.15  and  0.30  showed  similar  tendencies, 
n 

From  the  figures,  it  is  clear  that  there  are  some  differences.  Figure  5 
which  shows  the  added  mass  coupling  coefficient  A^  of  sway  into  yaw  also 
has  included  results  of  experiments  conducted  by  Vugts  (1970) . The  experi- 
ments were  forced  motion  tests  on  a 10  foot  segmented  model.  The  sectional 
added  mass  components  were  added  to  give  the  results  shown  in  the  figure.  The 


comparison  between  theory  and  experiment,  at  least  for  the  coefficient  of  A^  > 
;<'oms  reasonable.  (Recall  that  the  expressions  A^  and  B^.  ^ are  unaffected 
by  the  antisymmetry  and  constant  mid  body  assumptions.  Consequently,  they 


should  be  more  applicable  to  ship  shapes  than  any  of  the  coefficients  resul- 


ting from  a coupling  with  roll.) 

In  Figures  9 and  10,  the  ratios  of  t(?^  to  are  plotted  as 


jj.  «j 


, -«■ mm  %*  t; 
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functions  of  io/L/g  . These  ratios  are  speed  independent  since  both  T^ 

jk 

and  t!')  vary  linearly  with  velocity.  From  Table  1,  it  is  clear  that  the 
3 K 

real  part  of  corresponds  to  the  velocity  dependent  part  of  the  Salvesen, 

Tuck  and  Faltinsen  added  mass  and  the  imaginary  part  of  T^  corresponds 

to  the  velocity  dependent  part  of  their  damping  coefficient.  Since  T f 2 ^ 

jk 

represents  the  free  surface  integrals,  we  can  see  from  these  plots  the  rela- 
tive importance  of  these  terms  to  the  usual  forward  speed  terms  as  given  by 
Salvesen,  Tuck,  and  Faltinsen  (1970).  For  some  frequency  ranges,  t(^  , 

the  free  surface  integral,  is  equal  to  or  larger  than  the  usual  forward  speed 


theory 


Added  mass  coefficients  A and  A Figure  8 Damping  cofficients  B and  B 


APPENDIX 


FORCE  AND  MOMENT  ON  THE  HULL 

Since  most  of  the  details  for  finding  the  force  and  moment  on  the  hull 
due  to  horizontal  motions  are  similar  to  those  for  vertical  motions  as  presented 
by  Ogilvie  and  Tuck  (1969) , only  a brief  summary  of  the  derivation  will  be 
given  here. 

The  time  dependent  solution  is  given  by  the  ij;(x,y,z,t)  potential.  It 
must  satisfy  equations  (5)  - (7)  which  are  repeated  here  for  convenience 


ib  + \b  =0  in  the  fluid  domain 
yy  zz 

-oi2ip  + gip  = -iu)U(2ij/  + 2y  >b  + x ) on  z=0 
z x y y yy 


(5) 

(6) 


and  on  the  hull 


f£  = iam2(n2+  xn6)  + i^Cyn^  zn2)  + Un^ 

-U(n2+  xn(.)(n2Xzz+  n3xyz)  + Un4[(nVxv- 


3Ay  2 z 


+ Z(n2xyy+  n3Xyz) ~ y(n2Xyz+  n3Xzz)  ^ 


(7) 


There  is  also  a radiation  condition  that  must  be  satisfied.  More  will  be  said 
about  this  later. 


In  order  to  simplify  the  solution  of  the  forced  oscillation  potentials, 
Ogilvie  and  Tuck  (1969)  defined  the  following  quantities: 

n.  and  m.  for  j=2,4,6  where 

3 3 J 

n2=  n2 

n4=  yn3-  zn2 

V Xn2 


m2=  -n2Xyy-  n3Xyz 


V xyn3~  XzV  Z(n2xyy+  "aV  " y(n2Xyz+  n3Xyz} 


m6m  n2~  Xm2 
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Let  the  mean  hull  surface  be  defined  by  the  equation 


SQ(x,y,z)  = 0. 


Then  define  a potential  where 


'’j  + ‘•‘j  = 0 in  the  fluid 

Jyy  Jzz 


3<f> . 


= n . on  S - 0 
3N  j 0 


and 


— oj?4>  • + g<f>  . =0  on  z=0  . 

-*  Jz 


Note  that  = x <t>2 

Next  define  a potential  I*  where 

1* . . V . =0  in  the  fluid 

1 + 3 

yy  Jzz 


3 V. 


-r-~-  = m.  on  S.  = 0 

3N  3 0 


and 


+ gf . =0  on  z=0 

3 3z 


Finally  define  a potential  where 


£2 . + fi . =0  in  the  fluid 

3yy  3zz 


an. 

3N 


0 on  S0  = 0 


and 


+ gflj  = ~ (2d>  ^ + 2xy4>j  + Xyy^j)  on  z=0 

z x y 


r U 


Then  il>=r.  ( iu><p  + uy.-  a)2on.)n.  satisfies  all  the  equations  in  the  boundary 
j J J J J 

value  problems  (equations  (5)  - (7))  to  the  order  considered.  The  radiation 
conditions  for  the  4> ^ and  problems  are  the  usual  ones  which  require 

out  going  waves.  The  n..  potential  represents  an  anti-synmetric  pressure 
distribution  on  the  free  surface,  and  as  shown  by  Troesch  (1975),  its  radiation 
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condition  is  represented  by  linearly  growing,  anti-symmetric  out-going  waves. 

As  shown  by  Ogilvie  and  Tuck  (1969) , the  pressure  can  be  found  from 
Bernoulli's  equation  and  then  integrated  over  the  mean  position  of  the  ship's 
hull  to  give  a total  force  consistent  to  an  order  of  0(e3//26)  . It  has  the 
following  form: 

F . (t)  = -p  ds{>]  |"n  . (iw)  2 - Um . (iw)~j  <|>  n.+  Un  . E | (iio)  f 

J J k ~ 3 J K K j — k 


+ (iw)  3f2k  J nk> 

As  done  in  the  main  text,  let  us  define  several  quantities  as  follows: 


pj(t)  = +T]k’>k(t)  • 

T<°>  = -p(iw)2  ds  n.<f>, 

:k  : k 

so 

Tjk)  = -piwU  ds  (n  jH'k^j^k) 
s0 

t(2)  = -p (iw) 3U  ds  n.fi,  . 

]k  3 k 


See  the  main  text  for  the  complete  expressions  for  T_^  when  j=2,4,6  and 
k=2,4,6  . 

Consider  now  . The  far  field  behavior  of  <j>.  and  ¥.  are  given 

J J 

as 

, ..  , . vz  i(wt  - v I y I ) 
iwcfijOj  -*■  sgn(y)Aj(x)e  e 11 


sgn^Bj  (*)«VZe*  V y ^ 

as  y -+  oo  . Here  v is  the  wave  number  given  by  w2=vg  . Using  this  fact 
and  Green's  theorem  for  two  dimensions,  we  can  follow  the  method  used  by 
Ogilvie  and  Tuck  (1969)  to  easily  show  that 

t(i)  = 0 
T j j 


r 

r.  *'!  ■ 
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Tjk  = -,)iu,U  Jds(<(,jmk  " W 


The  following  expressions  are  given  for  specific  modes  of  motion: 


,u)  = _To)=  _pia)U  Jds  n^2 
S0 

24)  = _T42)  = "piwU  jds((t,2m4  ~ *4*2* 


T46)  ~ “T64)  = _pi(ijUL  jdsx(m2<()4  - m^)  + ds  n2<f>4 


To  simplify  T?j^  / we  need  the  far  field  behavior  of  the  ft  potential. 
Using  the  method  shown  by  Troesch  (1975)  for  a anti-symmetric  pressure  dis- 
tribution on  the  free  surface  we  can  write 

-«2UOj.y  9,n(y)e',z+i<“t-vlyl>[cj(x)  - A,  (x)(2-i|y|)] 


as  y 


Then  we  can  show,  using  a method  similar  to  that  used  by  Ogilvie  and  Tuck 
(1969) , that 

t(?)  =o 
jj 


T^|*  - - Pile.)3  ~ jdx["  | dy(<,|  - A 


l - A|e-2iv'') 


L 1, 


- h >¥ 


and  for  roll,  when  either  j=4  or  k=4  and  the  other  subscript  equals  2 or  6 
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-Tkj}  = PUW)3  f ldXfjdy[(<))j\  " Aj\xe‘2lVy) 

L yo  " 

' AjxV'2lVY)  + xy(Vk/  ♦jy+k)] 

* & e'2iuv°(AjX\  - aax] 


where  the  subscript  of  x means  the  derivative  of  the  subscripted  variable 
with  respect  to  x 

Note  the  following: 


i)  For  symmetric  hull  forms  m2 
with  respect  to  x 


m.  and  v are  odd  functions 
4 Ay 


ii) 


In  regions  of  constant  cross  sections, 
A ^ are  equal  to  zero. 


X , and 

y Jx 


iii)  In  the  end  sections  the  beam  of  the  hull  becomes  higher  order 

than  the  beam  at  mid  ships.  As  a result  the  oscillation  poten- 
tials <J>2  and  <t>4  , also  are  higher  order  there  than  in  the 

mid  sections. 


NUMERICAL  METHOD  FOR  DETERMINING  THE  ADDED  MASS  AND  DAMPING  COEFFICIENTS 

The  added  mass  and  damping  coefficients  were  found  in  two  steps.  First, 
the  potential  valid  on  the  hull  surface  was  determined.  Next,  a multipole 
expansion,  valid  outside  some  radius  R which  inclosed  the  body,  was  matched 
to  the  first  potential.  Using  tnis  scheme  significantly  lowered  the  computing 
time  of  evaluating  the  free  surface  integrals. 

The  potential  in  the  near  field  was  given  by  an  integral  representation 
of  sources  distributed  over  the  hull  surface.  The  source  strength,  Oj(C,n)  , 
was  found  by  solving  the  following  integral  equation: 

3*j(y'z)  ( 9G 

= ~1t0 j (y,z)  + jdfco^u.n)  (y,z,K#n) 

CH 

on  the  hull,  where  is  the  normal  velocity  on  the  hull  and  G(y,z,;C#n) 

is  a Green's  function  given  by 
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G(y,z;<;,n)  = log/(y-r,)?  + (z-n)?  - log/(y-02  + (z+n)2 
-2eV(z+n)Rfi(eiv|y’c|E1[v(z+n)  + iv | y-C | ]) 


A 


+27TieV(z+n)  - iVly_r'l 

and  E^(x+iy)  is  the  complex  exponential  integral.  The  method  for  solving 
the  integral  equation  is  given  in  Troesch  (1975) . 

Once  the  source  distribution,  o.(C,n)  was  known,  a circle  of  radius  R, 
where  R enclosed  the  entire  station  being  considered,  was  found.  Then  the 
potential  given  by  the  source  distribution  on  the  hull  was  matched  to  a multi- 
pole expansion  consisting  of  a dipole  and  wave  free  potentials  that  were  odd 
with  respect  to  the  x-z  plane.  The  matching  took  place  on  the  circle.  The 
multi-pole  expansion  had  the  following  form: 

<f>j(y.z)  = A j |~-  ^ eVZsgn(y)Re  |ielv  ^ y ^ (vz+iv  | y | ) J 


l y 

+ , + sgn(y)e 

n\> 


vz-iv|y| ) 


N-l 


2m+l f sin(2m-H) 9 . sin2m8 


+ * Ajm(vR)  "(vr)  2m+l  + Xj^pml 

m=i 

where  A_.  and  A are  coefficients  determined  from  the  matching  processes, 
r equals  /y2+z2  , and  0 is  the  angle  between  r and  the  negative  z 
axis. 

The  free  surface  integrals,  defined  as 


dy (4>|  - A2e"2lvy  ) 


*0 


and 


r 


dy(^^  - V4e’2ivV) 


2 4 


have  integrands  that  oscillate  with  a period  of  n . The  subroutine  that 
evaluated  the  integrals  used  Simpson's  rule  on  25  points  for  each  interval  of  tt . 
The  integration  was  terminated  when  a given  interval  made  no  significant  con- 
tribution to  the  total  integral. 
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& Engineering  Station  (Library) 

1 

Norfolk  Naval  Shipyard 

Newport,  RI  02840 

Technical  Library 

) 

Portsmouth,  VA  23709 

i 

Commanding  Officer  (L31) 

Naval  Civil  Engineering  Laboratory 

1 

Philadelphia  Naval  Shipyard 

Port  Hueneme,  CA  93043 

Philadelphia,  PA  19112 

Attn:  Code  240 

3 

Commander 

Naval  Ocean  Systems  Center 

1 

Portsmouth  Naval  Shipyard 

San  Diego,  CA  92152 

Technical  Library 

Attn:  Dr.  A.  Fabula  (4007) 

Portsmouth,  NH  03801 

Dr.  J.  Hoyt  (2501) 

Library  (13111) 

1 

Puget  Sound  Naval  Shipyard 

Engineering  Library 

, 

1 

Director 

Bermerton,  WA  98314 

Naval  Research  Laboratory 

Underwater  Sound  Reference  Division 

1 

Long  Beach  Naval  Shipyard 

P.0.  Box  8337 

Technical  Library  (246L) 

Orlando,  FL  32806 

Long  Beach,  CA  90801 

1 

Library 

1 

Mare  Island  Naval  Shipyard 

Naval  Underwater  Systems  Center 

Shipyard  Technical  Library 

Newport,  RI  02840 

Code  202.3 

Vallejo,  CA  94592 

1 

1 

Research  Center  Library 

Waterways  Experiment  Station 

1 

Assistant  Chief  Design  Engineer 

Corps  of  Engineers 

for  Naval  Architecture  (Code  250) 

P.0.  Box  631 

Mare  Island  Naval  Shipyard 

Vicksburg,  MS  39180 

Vallejo,  CA  94592 

2 

National  Bureau  of  Standards 

2 

U.S.  Naval  Academy 

Washington,  DC  20234 

Annapolis,  MD  21402 

Attn:  P.  Klebanof f (FM  105) 

Attn:  Technical  Library 

Fluid  Mechanics 

Dr.  Bruce  Johnson 

Hydraulic  Section 

3 

Naval  Postgraduate  School 

1 

AFFOL/FYS  (J.  Olsen) 

Monterey,  CA  93940 

Wright  Patterson  AFB 

Attn:  Library,  Code  2124 

Dayton,  OH  45433 

Dr.  T.  Sarpkaya 

Prof.  J.  Miller 

1 

Dept,  of  Transportation 

Library  TAD-491.1 

1 

Study  Center 

400  - 7th  Street,  S.W. 

National  Maritime  Research  Center 

Washington,  DC  20590 

U.S.  Merchant  Marine  Academy 

* 

S I 

Kings  Point,  L.I.,  NY  11024 

1 

T 
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1 U.S.  Merchant  Marine  Academy  2 

Kings  Point,  L.I.,  NY  11024 
Attn:  Academy  Library 

1 Library 

The  Pennsylvania  State  University 
Applied  Research  Laboratory  1 

P.0.  Box  30 

State  College,  PA  16801 

1 Bolt,  Beranek  & Newman 

50  Moulton  Street 

Cambridge,  MA  02138  1 

Attn:  Library 

1 Bethlehem  Steel  Corporation 

Center  Technical  Division 
Sparrows  Point  Yard 

Sparrows  Point,  MD  21219  1 

1 Bethlehem  Steel  Corporation 

25  Broadway 
New  York,  NY  10004 

Attn:  Library  (Shipbuilding)  1 

1 Cambridge  Acoustical  Associates,  Inc. 

1033  Mass  Avenue 

Cambridge,  MA  02138  1 

Attn:  Dr.  M.  Junger 

1 CALSPAN  Corporation 

P.O.  Box  235 
Buffalo,  NY  14221 

Attn:  Dr.  A.  Ritter  1 

Aerodynamics  Res.  Dept. 

1 Mr.  V.  Boatwright,  Jr. 

R&D  Manager  1 

Electric  Boat  Division 
General  Dynamics  Corporation 
Groton,  CT  06340 

1 Gibbs  & Cox,  Inc.  1 

21  West  Street 
New  York,  NY  10006 
Attn:  Technical  Info.  Control 

2 

1 Hydronautics,  Inc. 

Pindell  School  Road 
Howard  County 
Laurel,  MD  20810 
Attn:  Library 


McDonnell  Douglas  Aircraft  Co. 

3855  Lakewood  Blvd 
Long  Beach,  CA  90801 
Attn:  J.  Hess 

T.  Cebeci 

Lockheed  Missiles  & Space  Co. 

P.O.  Box  504 
Sunnyvale,  CA  94088 
Attn:  Mr.  R.L.  Waid,  Dept.  57-74 
Bldg.  150,  Facility  1 

Newport  News  Shipbuilding  & 

Dry  Dock  Company 
4101  Washington  Avenue 
Newport  News,  VA  23607 
Attn:  Technical  Library  Dept. 

Nielsen  Engineering  & Research  Inc.  ‘ 

510  Clude  Avenue 
Mountain  View,  CA  94043 
Attn:  Mr.  S.  Spangler 

J 

Oceanics,  Inc. 

Technical  Industrial  Park 
Plainview,  L.I.,  NY  11803 

Society  of  Naval  Architects 
and  Marine  Engineers 
One  World  Trade  Center,  Suite  1369 
New  York,  NY  10048 
Attn:  Technical  Library 

Sun  Shipbuilding  & Dry  Dock  Co. 

Chester,  PA  19000 

Attn:  Chief  Naval  Architect 

Sperry  Systems  Management  Division 
Sperry  Rand  Corporation 
Great  Neck,  NY  11020 
Attn:  Technical  Library 

Stanford  Research  Institute 
Menlo  Park,  CA  94025 
Attn:  Library  G-021 

Southwest  Research  Institute 
P.O.  Drawer  28510 
San  Antonio,  TX  78284 
Attn:  Applied  Mechanics  Review 
Dr.  H.  Abramson 
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1 Tracor , Inc . 

6500  Tracor  Lane 
Austin,  TX  78721 

1 Mr.  Robert  Taggart 

9411  Lee  Hgwy,  Suite  P 
Fairfax,  VA  22031 

1 Ocean  Engr  Department 

Woods  Hole  Oceanographic  Inst. 
Woods  Hole,  MA  02543 


l Colorado  State  University 

Foothills  Campus 
Fort  Collins,  CO  80521 
Attn:  Reading  Room, 

Engr  Res  Center 

1 University  of  California  at 

San  Diego 

La  Jolla,  CA  92038 
Attn:  Dr.  A . T.  Ellis 

Dept,  of  Applied  Math 


1 Worcester  Polytechnic  Inst. 

Alden  Research  Laboratories 
Worcester,  MA  01609 
Attn:  Technical  Library 

1 Applied  Physics  Laboratory 

University  of  Washington 
1013  N.E.  40th  Street 
Seattle,  WA  98105 
Attn:  Technical  Library 

1 University  of  Bridgeport 

Bridgeport,  CT  06602 
Attn:  Dr.  E.  Uram 

1 Cornell  University 

Graduate  School  of  Aerospace  Engr 
Ithaca,  New  York  14850 
Attn:  Prof.  W.R.  Sears 

4 University  of  California 

Naval  Architecture  Department 
Colleqe  of  Engineering 
Berkeley,  CA  94720 
Attn:  Library 

Prof.  W.  Webster 
Prof.  J.  Paul ling 
Prof.  J.  Wehausen 


2 Florida  Atlantic  University 

Ocean  Engineering  Department 
Boca  Raton,  FL  33432 
Attn:  Technical  Library 

Dr.  S.  Dunne 

2 Harvard  University 

Pierce  Hall 
Cambridge,  MA  02138 
Attn:  Prof.  G.  Carrier 

Gordon  McKay  Library 

1 University  of  Hawaii 

Department  of  Ocean  Engineering 
2565  The  Mall 
Honolulu,  HI  96822 
Attn:  Library 

1 University  of  Illinois 

Urbana,  IL  61801 
Attn:  Dr.  J.  Robertson 

3 Institute  of  Hydraulic  Research 

The  University  of  Iowa 

Iowa  City,  IA  52240 
Attn:  Library 

Dr.  L.  Landweber 
Dr.  J.  Kennedy 


3 California  Institute  of  Technology  1 

Pasadena,  CA  91109 
Attn:  Aeronautics  Library 
Dr.  T.Y.  Wu 
Dr.  A.J.  Acosta 

1 

1 Docs/Repts/Trans  Section 

Scripps  Institution  of 
Oceanography  Library 
University  of  Calif.,  San  Diego 
P.0.  Box  2367 
La  Jolla,  CA  92037 


The  John  Hopkins  University 
Baltimore,  MD  21218 
Attn:  Prof.  0.  Phillips 

Mechanics  Dept 

Kansas  State  University 
Engineering  Experiment  Station 
Seaton  Hall 
Manhattan,  KS  66502 
Attn:  Prof.  D.  Nesmith 


T 
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1  University  of  Kansas 

Chm  Civil  Engr  Dept  Library 
Lawrence,  KS  60644 

1 Fritz  Engr  Laboratory  Library 

Depart  of  Civil  Engr 
Lehigh  University 
Bethlehem,  PA  18015 

3 Department  of  Ocean  Engineering 

Massachusetts  Institute  of 
Technology 
Cambridge,  MA  02139 
Attn:  Prof.  P.  Mandel 

Prof.  M.  Abkowitz 
Dr.  J.  Newman 

1 Engineering  Technical  Reports 
Room  10-500 

Massachusetts  Institute  of 
Technology 

Cambridge,  MA  02139 

5 St.  Anthony  Falls  Hydraulic 

Laboratory 

University  of  Minnesota 
Mississippi  River  at  3rd  Avenue  SE 
Minneapolis,  MN  55414 
Attn:  Dr.  Roger  Arndt 

Mr.  J.  Wetzel 
Prof.  E.  Silberman 
Mr.  J.  Killen 
Dr . C . Song 

3 Department  of  Naval  Architecture 

and  Marine  Engineering 
North  Campus 

The  University  of  Michigan 
Ann  Arbor,  MI  48109 
Attn:  Library 

Dr.  T.F.  Ogilvie 
Prof.  F.  Hairanitt 

2 College  of  Engineering 
University  of  Notre  Dame 
Notre  Dame,  IN  46556 
Attn:  Engineering  Library 

Dr.  A.  Strandhagen 

2 New  York  University 

Courant  Inst,  of  Math.  Sciences 
251  Mercier  Street 
New  York,  NY  10012 
Attn:  Prof.  A Peters 
Prof.  J.  Stoker 


1 Department  of  Aerospace  & 

Mechanical  Sciences 
Princeton  University 
Princeton,  NJ  08540 
Attn:  Prof.  G.  Mellor 

3 Davidson  Laboratory 

Stevens  Institute  of  Technology 
711  Hudson  Street 
Hoboken,  NJ  07030 
Attn:  Library 

Mr.  J.  Breslin 
Mr.  S.  Tsakonas 

1 Applied  Research  Laboratory 

Library 

University  of  Texas 
P.0.  Box  8029 
Austin,  TX  78712 

1 College  of  Engineering 
Utah  State  University 
Logan,  UT  84321 
Attn:  Dr.  R.  Jeppson 

2 Stanford  University 
Stanford,  CA  94305 

Attn:  Engineering  Library 

Dr.  R.  Street 

3 Webb  Institute  of  Naval  Architecture 

Crescent  Beach  Road 

Glen  Cove,  L.I.,  NY  11542 
Attn:  Library 

Prof.  E.V.  Lewis 
Prof.  L.W.  Ward 

1 National  Science  Foundation 

Engineering  Division  Library 
1800  G.  Street  N.W. 

Washington,  DC  20550 

1 University  of  Connecticut 

Box  U-37 

Storrs,  CT  06268 
Attn:  Dr.  V.  Scottron 

Hydraulic  Research  Lab 

1 Applied  Research  Laboratory 

P.0.  Box  30 

State  College,  PA  16801 
Attn:  DR.  B.  Parkin,  Director 
Garfield  Thomas  Water 
Tunnel 
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Dr.  Michael  E.  McCormick 
Naval  Systems  Engineering  Department 
U.S.  Naval  Academy 
Annapolis,  MD  21402 

Dr.  Douglas  E.  Humphreys  (Code  794) 
Naval  Coastal  Systems  Laboratory 
Panama  City,  FL  32401 

Dr.  Allen  H.  Magnuson 
Aerospace  and  Ocean  Engineering  Dept. 
Virginia  Polytechnic  Institute  and 
State  College 
Blacksburg,  VA  24061 
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